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Abstract. Magnetic helicity effects are discussed in laboratory and astrophysical 
settings. First, dynamo action in Taylor-Green flows is discussed for different boundary 
conditions. However, because of the lack of scale separation with respect to the 
container, no large scale field is being produced and there is no resistively slow 
saturation phase as otherwise expected. Second, the build-up of a large scale field is 
demonstrated in a simulation where a localized magnetic eddy produces field on a larger 
scale if the eddy possesses a swirl. Such a set-up might be realizable experimentally 
through coils. Finally, new emerging issues regarding the connection between magnetic 
helicity and the solar dynamo are discussed. It is demonstrated that dynamos with a 
non-local (Babcock-Leighton type) alpha effect can also be catastrophically quenched, 
unless there are magnetic helicity fluxes. 



PACS numbers: 52.30.Cv, 47.65.Md, 96.60.Hv 
1. Introduction 

Many astrophysical dynamos are driven by helical flows via an a effect. Helical flows 
are also employed in all laboratory realizations of liquid metal dynamos. One of the 
remarkable properties of helical dynamos is that they produce large scale fields that are 
helical. Since net helicity is conserved, this can only happen if here is a production of an 
equal amount of small scale fields with opposite helicity This has been demonstrated 
through various numerical studies (Brandenburg 2001, Mininni et al. 2005a). 

An obvious question is whether these properties can also be seen in experimental 
realization of fluid dynamos. The preliminary answer to this is no, because there is no 
scale separation. We will return to this in more detail and present new calculations of 
Taylor-Green flows that model the flow in the French VKS2 experiment (Monchaux et 
al. 2007). To address the issue of the lack of scale separation we also study a model 
designed to show the development of large scale magnetic fields that are driven from a 
small localized source. Finally, we discuss some new issues regarding the solar dynamo. 

We assume some level of familiarity with the concept of magnetic helicity 
conservation and the resistively limited saturation phenomenon found in the nonlinear 
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evolution of large scale dynamos. In this connection we highlight the usage of the word 
"catastrophic" , which is meant to indicate that the result depends on the value of the 
magnetic Reynolds number. This applies in particular to the a effect which is now 
known to be catastrophically quenched when there is no flux of magnetic helicity. A 
recent review on the development of the last five years can be found in Brandenburg & 
Subramanian (2005a). 

Throughout this paper we present original results that have not yet been presented 
earlier. Many of the cases considered are motivated by the recent developments in 
laboratory dynamos and solar dynamo modeling. Many of the simulations have been 
carried out at relatively low resolution. The present results are therefore tentative 
and need to be followed up using higher resolution simulations. However, the results 
presented here do reflect the current state of affairs in this field, which is indeed one of 
the objectives of this paper. 



2. Governing equations 

We consider here the isothermal and weakly compressible case, which means that the 
pressure gradient term can be written as a gradient of the pseudo enthalpyjij] h — In p, 
where c s = const is the isothermal speed of sound and p is the density. The induction 
equation can then be written in an analogous form in terms of the magnetic vector 
potential by using the pseudo Lorenz gauge with a freely specified speed c^, which is 
for practical reasons less than the speed of light, and may well be equal to the speed of 
sound c s . The set of equations is thus 

dA dd> 

— = -V<p + £ cxt - V J + uxB, _^ = _ c Jv-A, (1) 

^ = -Vh + F ext -vQ + Jx B/p, ^ = -c s 2 V • U, (2) 

where U is the velocity, B = V x A is the magnetic field, expressed in terms of the 
magnetic vector potential A, J = Vx VxA is the current density, Q = V x V x U 
is the double curl of the velocity, D/Dt = d/dt + U ■ V is the advective derivative, and 
£ cxt and F ex t are external forcing functions, to be specified later. The current density 
is measured in units where the vacuum permeability is unity. The main analogy we 
want to stress here is that between the electric potential in the uncurled induction 
equation and the pseudo enthalpy h in the momentum equation. We should point out 
that, in order to not disturb the analogy between the two equations, we have ignored a 
correction term in equation ([1]) where Q has to be replaced by 

Q^Q + VV-C/ + S- Vln(pi/), (3) 

where Sy = \{Uij + Ujj) — ■ U is the traceless rate of strain tensor. For 

incompressible flows these correction terms vanish, and they are small for weakly 

% For a polytropic gas the enthalpy can be written as h = (7 — l)^ 1 c^ (p/ po) 7_1 > which reduces to 
Cg \n(p/po) + const in the limit of 7 —> 1. In the following we include the above constant in our 
definition of the enthalpy and refer to it therefore as the pseudo enthalpy. 
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compressible (small Mach number) flows considered here. Nevertheless, in the numerical 
computations these extra terms in equation ([3]) are always included. 

If we think of the velocity field being primarily just the vector potential for the 



vorticity (see Appendix A), then we can regard the continuity equation, Dh/Dt = 
— CgV • U, as the gauge condition for the velocity in the Lorenz-like gauge. In this 
sense we can interpret the Lorenz gauge condition for A as a continuity equation for <fi. 
Note, however, that A is not a physically measurable quantity. 

3. Lorenz versus Weyl gauge 

Magnetic and electric fields are invariant under the gauge transformation 

A' = A + VA, (4) 

For numerical purposes it is often convenient to choose the gauge A = f <f)dt, which 
implies that <// = 0. Thus, instead of equation (jTJ one solves the equation dA' /dt = —E. 
The latter is usually referred to as the Weyl gauge. 

In this section we compare the two gauges. Obviously the magnetic field is the 
same in both cases, but the magnetic vector potential is not. In Fig. [I] we compare the 
results obtained with the two techniques in the case of a dynamo driven by fully helical 
turbulence. The detailed forcing function used here is given in §Appendix B We choose 



32 3 mesh points, a forcing wavenumber of k{ = 3k%, where ki is the smallest wavenumber 
in the box, and a forcing amplitude of 0.07, which results in an rms velocity of about 0.2. 
Viscosity and magnetic diffusivity are chosen to be 5 x 10~ 3 , so the magnetic Reynolds 
number is R m = u rms / (rjkf) = 13. For details about these and similar runs at larger 
magnetic Reynolds number, see Brandenburg (2001), which is also where the resistively 
limited saturation phase with 

(B 2 ) oc 1 - e" 2 " fc i (i "* sat) (6) 

was proposed. Here, £ sat is the time where the small scale magnetic field saturates. The 
dotted line in the lower left panel shows this behavior. 

For the present studies we usedjust 32 3 meshpoints, which can easily run on one 
processor using the Pencil C0DE3I, which is a non-conservative, high-order, finite- 
difference code (sixth order in space and third order in time) for solving the compressible 
hydrodynamic equations. The particular run presented in this section is actually one of 
the sample runs (helical-MHDturb) that come with the code. Here and in all other cases 
presented in this paper we use c s = = 1. 

Obviously, the resistive saturation phase becomes more prominent at larger 
Reynolds numbers. However, the point of this simulation was to illustrate the differences 
between Lorenz and Weyl gauges. Qualitatively, it appears that the magnetic vector 

§ http : //www . nordita. dk/sof tware/pencil-code 
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Figure 1. Comparison of runs in Weyl and Lorenz gauges. In this simulation a large 
scale magnetic field develops owing to the helicity of the velocity forcing function. 
While the magnetic field is the same in both gauges (second and third column in the 
second row) and vector potential is different (second and third column in the first row). 
These slices are taken after 3 magnetic diffusion times. For comparison, the evolution 
of magnetic energy in the total field (upper left panel) and in the mean field (lower 
left panel) are shown, where time is given in diffusion times. 



potential in the Lorenz gauge is smoother than that in the Weyl gauge. In general, 
the divergence of the magnetic vector potential is small, so for all practical purposes 
the Lorenz gauge is close to the Coulomb gauge. One may hope that in the presence of 
open boundary conditions, where a gauge-invariant magnetic helicity is harder to define, 
the Lorenz-gauged magnetic vector potential may provide some meaningful guidance 
regarding the escape of magnetic helicity density and the associated magnetic helicity 
fluxes (cf. Subramanian & Brandenburg 2006). 

4. Taylor-Green flow dynamos 

The Taylor-Green (TG) flow is often studied in connection with the von Karman Sodium 
(VKS) dynamo experiment in Cadarache in France (Monchaux et al. 2007). The TG 




Figure 2. Cross-section through the middle of the domain (y — it/2 in the upper row 
and z = ir/2 in the lower row) of velocity (left hand column) and magnetic field (right 
hand column). Velocity and magnetic fields are represented by vectors in the plane 
and colors (gray scales) indicate the normal components. The color (gray scale) bar 
indicates the contour range symmetrically about zero within the maximum values of 
either sign. The light red color (intermediate gray) corresponds to zero. All boundaries 
are assumed to be perfectly conducting. Re=5 and R m = 1000. 



cos k x sin k y cos k z I, — — <x,y, z< — . (7) 

k k 



flow is given by u = J- ext / (ukf ), where k{ = V3k and 

( + sin k Q x cos koy cos k§z 1 

J~ cxt — 2/o 

V o 

The forcing function is normalized such that (^"ext) 1 ^ 2 = fo- There is a vast literature on 
this flow. Due to a large number of symmetries allowing computational simplifications, 
large Reynolds numbers can be achieved. This flow has therefore been used to study 
singularities and turbulence (Nakano 1985, Brachet 1991). In more recent years dynamo 
action for this flow has been studied (Nore et al. 1997), especially at low magnetic 
Prandtl numbers (Ponty et al. 2004, 2005, Mininni et al. 2005b). 

In the following we use fo = 0.006, v = 0.02, so that M rms = 0.1. We adopt units 
of length where k$ = 1. In agreement with common practice, we use in this section 




Figure 3. Same as Fig. [3j but with vertical field boundary conditions at z = and 
z = 7T. Re=5 and E m = 200. 



Reynolds numbers based on the scale l/k rather than 1/kf, i.e. Re = u rms /(vki) 
and R m = u rms /(r]ki). No magnetic forcing is applied, i.e. £ ext = 0. The VKS flow 
corresponds to one eight of the full domain, i.e. < x, y, z < tt. Since dynamo action in 
the TG flow is normally studied in triply periodic domains we show in Fig. [2] numerically 
obtained results for a domain bounded by perfect conductors on all sides. In Fig. [3] we 
give the corresponding result for a domain where we adopt a 'vertical field' boundary 
condition on z = and z = tt. The latter is supposed to simulate the experimental 
situation of soft iron boundaries with a large permeability (Fauve & Petrelis 2003, 
Morin 2005, E. Dormy, private communication; see also Kenjerescaron & Hanjalicacute 
2007). Finally, in Fig. H] we show the result for a triply periodic domain. The resulting 
saturation field strengths for the three cases are summarized in Table [TJ 

Evidently, the dynamo operation is quite sensitive to the boundary conditions in 
that the critical value of R m is about 5 times larger when all domain boundaries are 
perfectly conducting. In this light it is not too surprising that no dynamo was found 
before using the soft iron lids at the two ends of the VKS experiment (Monchaux et al. 
2007). Indeed, the model with the vertical field condition has the lowest critical magnetic 




Figure 4. Same as Fig. [3j but for triply periodic boundary conditions in the larger 
domain, — 7r < x,y,z < n. The box in the first quadrant marks the location of 
the computational domains used for the results shown in Figs [2] and [3] Note that 
the magnetic field is concentrated preferentially along the cell boundaries. Re=5 and 
R m = 200. 



Table 1. Saturation values of B rms = [B 2 ) 1 / 2 for TG flow dynamos in triply-periodic 
domains of volume (27r) 3 (second column) compared with those in one eighth of the 
volume (7r 3 ), where the boundaries in the x and y directions are perfectly conducting, 
and those in the z direction either using a vertical field boundary condition (third 
column) or also perfectly conducting (last column). 

R m periodic vert, field perfect cond. 

200 0.027 

500 0.053 0.069 
1000 0.041 
2000 0.054 



Reynolds number of all three cases considered. The field shows a narrow vertical flux 
tube in the middle of the domain. Obviously, given that our fluid Reynolds number is 
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■q kf z t r/ kf z t 

Figure 5. Saturation behavior of a run with a vertical field boundary condition 
similar to that shown in Fig. [3] with Re = 5 and R m = 500. The dotted lines represent 
failed attempts to match the saturation behavior to the functional form of equation ([6]), 
so the saturation is not resistively limited in this case. Note the bursty oscillations in 
the electric potential difference during early times when the field becomes strong. 




77 kt z t 7] kf z t 

Figure 6. Solution of the one-mode truncation of the Blackman-Field model with 
R m = 500 and St = 20. The left hand panel more clearly shows the oscillations around 
the time rjkft = 0.35, while the right hand panel shows the full saturation phase, which 
is matched perfectly by equation ([6|) as shown by the dotted line. 

small, the velocity field is laminar. Larger Reynolds numbers could be achieved with 
more meshpoints. 

In Fig. Owe show that the root mean square magnetic field strength saturates after 
about one magnetic diffusion time. Nevertheless, there is no prolonged saturation phase 
as it was seen in Fig. [TJ or in other forced turbulence simulations where the forcing 
wavenumber kt/k\ is large. In Fig. [5] we also show the electric potential difference over a 
distance of 4 meshpoints. It turns out that the electric potential difference shows strong 
bursts during times when the field is strong. At late times the potential difference 
oscillates with a frequency of about 0.14, which is about 7 times smaller than the basic 
frequency, c^k\. 
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Figure 7. Comparison of the evolution of the rms values of the mean magnetic 
field obtained by averaging over xy planes (solid lines), yz planes (dotted lines), and 
xz planes (dashed lines) for domains of different size. Only in the case of a bigger 
domain, — 47r < x,y,z < 47r, is there a weak mean field (solid curve). In the case of 
the smaller domain, — n < x,y, z < w, there is in principle also a mean field, but it it 
noisy and there is no preferred plane of averaging. In both cases v = y = 2 x 1CP 4 . 



The nature of the oscillations in Fig. [5] remains unclear. Similar oscillations 
have been seen in solutions of mean-field dynamo equations derived under the r 
approximation where the explicit time derivative of the electromotive force is retained 
(Blackman & Field 2002). Such oscillations have also been seen in direct simulations 
of passive scalar turbulence, when the forcing scale is close to the scale of the domain 
(Brandenburg et al. 2004). This is actually the case here, because k^/ki = \/3 is 
close to unity. However, there appear to be two problems with this interpretation, 
which are illustrated using a numerical solution of an appropriately adjusted version 
of the Blackman-Field model; see Fig. [6j This model is here solved in the one-mode 
truncation; see §Appendix C Firstly, the frequency of the oscillations is in the model 
^osc = u vms k{/ 'y/3, but this is about 20 times higher than what is actually seen. Secondly, 
the slow saturation phase, as described by equation (jHJ), should still be seen in the 
saturation of the mean field. For these reasons we can conclude that the behavior seen 
in Taylor-Green flow dynamos is not described by standard mean field models. It should 
also be noted that any mean field in these dynamos cannot be seen as spatial averages. 

Temporal averages may be a useful way of analyzing the evolution of mean fields in 
Taylor-Green flow dynamos. It turns out that after saturation (rjkft > 3) the time 
averaged field contains about 88% of the total magnetic energy. At earlier times 
the oscillations are well reproduced by time averages over short enough time spans 
(rjkfAt = 0.04). No evidence for a slow resistively limited saturation phase is seen 
even for intermediate time averages. Again, this supports the idea that the magnetic 
fields seen in Taylor-Green flow dynamos cannot be described by the existing mean field 
approach. 
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In hindsight the absence of a slow resistively limited saturation phase is not too 
surprising because the wavenumber of the forcing function, kf = V3ko, is close to the 
box wavenumber, so there is no scale separation. This means that any large scale field, 
if it exists, could only be generated at almost the same scale. More importantly, because 
on short time scales no magnetic helicity can be generated, the positive magnetic helicity 
in the forcing must be balanced by the negative magnetic helicity at the same scale, 
leading essentially to a cancellation of the magnetic field. 

In order for there to be a large scale field, it is important to allow for a domain 
size that is sufficiently large to accommodate a field whose scale is at least a few times 
bigger than the eddy scale. In Fig. [7| we show the growth of the rms field strength of 
mean fields defined by averaging over different horizontal planes. Most interesting is 
the xy plane, and the corresponding averages are denoted by solid lines. Only in the 
case of the larger domain a mean field can be identified. In all other cases the resulting 
averages are dominated by 'noise'. Nevertheless, even in the case of the bigger domain 
the amplitude of the mean field is small compared with the equipartition field strength, 
B cq = (/i p^ 2 ) 1/2 . 

5. Inverse transfer from a localized source 

Given that the problem of not seeing a large scale field with a prolonged saturation phase 
of (B 2 ) is related to the lack of scale separation, we now devise a special experiment 
that can probably also easily be realized in the laboratory. The main idea is that the 
inverse transfer from small to large scales is primarily connected with the conservation 
of magnetic helicity (Frisch et al. 1975), so it is a magnetic phenomenon that is best 
demonstrated using magnetic forcing in the induction equation. This was done in the 
seminal paper by Pouquet et al. (1976) and the resistively slow build-up of magnetic 
energy has also been studied by Brandenburg et al. (2002; see their Fig. 12). 

The forcing needs to be helical in order to produce magnetic helicity. In order to 
have scale separation we apply the forcing only in small localized parts of the domain. 
Here we restrict ourselves to only a single localized source that has been modeled as 

£ ext (x) =(tV x (<pz) + V x V x (<pz), (8) 

where we have chosen a gaussian profile for <p with 

(p(x) = e exp(-x 2 /R 2 ) } (9) 

and a = 1 is chosen for a helical flow. This yields 

/ -2ay/R 2 + Axz/R A \ 
£ cxt (x)=\ 2ax/R 2 + Ayz/R 4 exp(-r 2 /i? 2 ). (10) 

\A(R 2 -x 2 -y 2 )/R A ) 
Localized forcing functions of this form yield a peak in the spectrum at 

fc peak = 2/R; (11) 
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Figure 8. Evolution of magnetic energy spectra and slices of the magnetic field, 
driven by a localized helical electromotive force at the center with a = 1 and 
eo = 0.02. The vertical arrows indicate the effective forcing wavenumber, as obtained 
from equation (fTTj) . Here, v = y = 5 x 10~ 3 , and so the three times shown in the figure 
correspond to ykft = 0.05, 2, and 3. 



see also Mee & Brandenburg (2006), where a potential momentum forcing proportional 
to Vy? was adopted. For all cases presented here we have chosen R = 0.5, so /c pea k = 4. 
The forcing amplitude is varied between eo = 0.005 and 0.02 such that the maximum 
flow and Alfven speeds remain subsonic. 

Looking at Fig. [8] we see that a magnetic eddy is produced that begins to swell up 
until it loses its original up-down orientation and tilts sideways in a somewhat irregular 
manner. At the same time the overall magnetic energy has decreased somewhat, but 
most of the spectral energy is now at large scales. This type of behavior is not seen 
without helicity injection, i.e. for a = 0, as is demonstrated in Fig. EO 



6. Connection with the solar dynamo problem 

We now wish to discuss the possible effects of magnetic helicity conservation that could 
be important for the solar dynamo. Obviously, only pieces of this question can be 
addressed at this point. Much of this has been discussed elsewhere (e.g. Blackman & 
Brandenburg 2002, 2003) and is summarized in various reviews (e.g. Brandenburg et 
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Figure 9. Magnetic energy spectrum and a slice of the magnetic field at t = 600 
for the case of a non-helical localized steady forcing at the center with a = and 
cq = 0.005. Again, v = i] = 5 x 10 -3 , and so the time shown in the figure corresponds 
to r\k\t = 3. 

al. 2002, Brandenburg & Subramanian 2005a). One would generally expect magnetic 
helicity conservation to be important for affecting the cycle amplitude unless magnetic 
helicity is allowed to escape the dynamo domain. 

Shear clearly plays an important role in dynamo processes. Two distinct effects can 
be identified. On the one hand shear can lead to so-called afl dynamo action which is 
often oscillatory. On the other hand, shear can allow local magnetic and current helicity 
fluxes along lines of constant angular velocity (Vishniac & Cho 2001, Subramanian & 
Brandenburg 2004, 2006). We focus here on the former effect. In a closed domain 
magnetic helicity conservation acts as to produce a "counterproductive" alpha effect 
that can saturate the dynamo, but it would not change the cycle frequency unless the 
turbulent magnetic diffusivity was also catastrophically affected. There has so far not 
been any clear evidence for catastrophic rjt quenching. Thus, we are at present not able 
to come to a conclusive results, which is mainly a result of limited computing resources 
available. However, it is also clear that in all cases the magnetic field is in strong 
excess of the kinetic energy, and so any quenching would be dominated by classical 
(non-catastrophic) effects. 

6.1. Catastrophic rj t quenching? 

For a saturated afl dynamo the cycle frequency is equal to uo cyc = r]^k\ (Blackman 
& Brandenburg 2002), where r^ T = r] t + 77 is the total (turbulent plus microscopic) 
magnetically quenched diffusivity and k\ is the largest wavenumber of the domain. This 
property is the key to what is perhaps the most robust method for determining the 
possible dependence of r] t on the magnetic field. 

This was already attempted in Brandenburg et al. (2002) using simulations with 



Magnetic helicity effects in astrophysical and laboratory dynamos 



13 




Figure 10. Visualization of the toroidal field B y on the periphery of the shearing 
box in the presence of forced turbulence for shear strengths S = —0.2, v = 5 x 10 -3 , 
and two different values of 77, 5 x 10 -3 on the left and 5 x 10~ 4 on the right. 

a sinusoidal shear profile, i.e. U y = U sin k\X. However, the result was not completely 
conclusive, because the solutions developed some additional complexity which resulted 
partly from the fact that shear strength and sign changed with x. In order to avoid this 
problem we consider now a linear shear profile of the form 

U = (0, Sx, 0) T , (12) 

where S is the shear rate. We solve the governing equations using a shearing box 
approach with shearing-periodic boundary conditions; see Hawley et al. (1995) and 
Brandenburg et al. (1995). 

Table 2. Summary of runs with uniform shear and different values of the magnetic 
diffusivity and either constant kinematic viscosity or constant magnetic Prandtl 
number. 
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The result of such simulations is shown in Fig. [10] for S = —0.2, v = 5 x 10~ 3 , and 
i] = 2 x 10~ 3 at a resolution of 64 3 (on the left) and 77 = 5 x 10~ 4 at a resolution of 256 3 
(on the right). In both cases the field is cyclic corresponding to an upward traveling 
wave, and the times are chosen such that the magnetic field is in approximately the 
same phase in both figures. The fact that the dynamo wave travels in the positive z 
direction is well understood as a consequence of a negative effective a (due to positive 
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Figure 11. Space-time, butterfly, or simply zt diagram of B x (z,t) and B y (z,t) for 
forced turbulence with shear, with S = —0.2, v — 5 x 1CP 3 and r\ — 5 x 10~ 3 . 



kinetic helicity in the forcing) and negative shear. The traveling wave behavior is best 
seen in a space-time diagram shown in Fig. [TTJ for the simulation with rj = 2 x 10~ 3 . 
It turns out that c<j cyc = 0.0032, and since ki = 1 this means 7/t = 0.0032 and hence 
?7t = 0.0012 or ?7 t /?7 = 0.6. This is still quite small, but we will be continuing to run 
models at larger resolution to increase the value of rjt/rj. In Table [2] we have summarized 
these parameters also for a few other runs. 



6.2. Catastrophic a quenching and non-locality 

The basic mechanism of catastrophic a quenching is now reasonably well understood. 
The basic recipe is this: whatever the mean turbulent electromotive force £ is, it leads 
not only to the production of large scale magnetic fields via the standard mean field 
equation, 

dB 

~cW 

but it also leads to the production of a magnetic alpha effect, «m, which characterizes 
the production of internal twist in the system. Its governing equation is 

^'^- 2 ^fi? + 3- (14) 



V x [U x B + £- rjiMoJ , (13) 
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Here we have allowed for the possibility of an additional flux of small scale magnetic 
helicity. Also, in this equation R m is meant to represent 77^/77 by definition; see Blackman 
& Brandenburg (2002) for details. Note that the magnetic helicity equation is unaffected 
by the large scale velocity term, U x B. 

It is important to realize that the possibility of catastrophic quenching is quite 
general and not restricted to the local a effect formula considered here. In some 
recent solar dynamo models the so-called Babcock-Leighton mechanism is used (e.g. 
Dikpati & Charbonneau 1999). Here one assumes that there is some source term to the 
electromotive force that is localized near the surface, but it is proportional to the toroidal 
field a the bottom of the convection zone. This effect is sometimes thought of being 
distinct from the usual a effect in that it allows for the generation of super-equipartition 
field strengths. 

Formally, the Babcock-Leighton mechanism is just a non-local a effect where the 
multiplication in aB is replaced by a convolution, ao B, where a is an integral kernel. 
Applying the same idea also to the turbulent magnetic diffusivity f)t(z, z'), leads to an 
expression for S(z,t) of the form 

£(z,t)= &{z,z')B(z',t)dz'- fj(z,z')J(z',t)dz'. (15) 

J Z\ J Z\ 

The possibility of nonlocal a and rj effects has been inferred also from simulation data 
of magneto-rotational turbulence in accretion discs (Brandenburg & Sokoloff 2002). 
Nonlocal a effects are conveniently characterized in terms of its spectral decomposition, 

rz2 

a(k) = / sin kz sin kz' a(z, z') dzdz'. (16) 



z\ 



This technique is particularly convenient because in Fourier space the convolution 
corresponds to a multiplication. 

In the following we restrict ourselves to a simple expression of the form 

a(z, z) = a g ou t{z) g m (z'), (17) 



where 



9in{z') = \ 



1 + erf 




(19) 



are simple profile functions representing the peak of the source function near z = z 2 with 
a sensitivity for fields located near z = Z\. For the following we choose — z\ = z 2 = 2.5 
and d = 0.05 in the domain —71 < z < n; see the right hand panel of Fig. [12j The 
resulting profile of a(k) according to equation fflBT) is shown in Fig. [13j 

We solve equations ( 1131) and ( |T^1) numerically for the case 7- = using an implicit 
scheme for «m by writing the equation in the form 

+ 2rf ^f + ^ \ + Rm WJ ^ + ^ = °' (20) 
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Figure 12. Profiles of the magnetic field and the magnetic contribution to the a 
effect for the nonlocal dynamo model with R m = 10 3 . 




Figure 13. Profile of the a kernel for the nonlocal (Babcock-Leighton type) dynamo 
model computed using equation (fTI)]) . 




200 400 600 800 1000 

Figure 14. Evolution of the magnetic energy of the large scale field for R m between 
10 2 and 10 4 for the nonlocal dynamo model. For large of i? m the linear growth 
rate reaches an asymptotic value, but the nonlinear saturation amplitude continues 
to depend on R m . 
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where = £ — ay[B is the electromotive force without the magnetic quenching term. 
Equation (120]) is then solved for at the new time, t n+ \ = t n + 5t, using for B the 
value at the present time level t n , i.e. 



1-Q) a$ L -2r h k$6t 



B 2 



(21) 



'eq 

where Q = i]kf5t(l + R m B 2 /B 2 q ). We have considered a model using linear shear of 
the form f[T2"j) . The strength of shear and a effect are quantified by the non-dimensional 
numbers 

S OL 

Cs = — rj, C a = — — , (22) 

where k% = 1 is the smallest wavenumber in the computational domain. In the following 
we use Cs = 100, C a = 0.1, and kf/ki = 5. In many cases an explicit treatment of the 
«m equation suffices (e.g. Blackman & Brandenburg 2002), but in the present case an 
explicit solution algorithm was found to be unstable. 

It turns out that most of the field is generated in the middle of the domain, 
while most of the quenching via «m occurs near the top layers around z = z 2 ; see 
Fig. [12j Nevertheless, as expected, this model experiences still catastrophic quenching; 
see Fig. [141 These results look quite similar to those obtained for local a profiles 
(Brandenburg & Subramanian 2005b). We note that for R m = 10 4 it is important to 
perform the calculations using double precision arithmetics. 



6. 3. Location of the dynamo 

The location of the solar dynamo is rather uncertain. Particularly unclear is the location 
where most of the toroidal field resides. The standard thinking since the 1980s is that 
most of the toroidal field can only reside at the base of the convection zone, because 
magnetic buoyancy would remove the field on a short time scale. However, it turned out 
that downward pumping of magnetic field helps to keep the magnetic field inside the 
convection zone. On the other hand, the toroidal field has been argued to be actually 
quite strong, such that it would exceed the equipartition field strength by factors as 
large as 100. Furthermore, stability of such strong fields even beneath the convection 
zone have recently been put under doubt (Arlt et al. 2007, Kitchatinov & Riidiger 2007). 
Whether or not such strong fields could be generated by a turbulent dynamo is unclear. 
This led to the proposal that the solar dynamo may instead be located in the bulk of 
the convection zone and that the mean fields are at most comparable in strength to the 
equipartition field strength (Brandenburg 2005). 

Astrophysical dynamos are not expected to be catastrophically quenched, because 
they are likely to shed excess small scale magnetic helicity through magnetic helicity 
fluxes. This would mean that the quenching through «m is much weaker than the 
other nonlinearities, e.g. suppression of the mean flow or of q;k itself. The latter can be 
approximated by the more conventional expression ok = «ko/(1 + B 2 /B 2 q ) (for detailed 
expressions see Rogachevskii & Kleeorin 2000). 
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If the solar dynamo does indeed work in a distributed fashion, then the meridional 
circulation is probably no longer important for determining the cycle period and the 
equatorward migration of the magnetic flux belts. Instead, the dynamo may be 
essentially of a standard aQ type. However, an important problem arises from the 
fact that the near-surface shear layer, where dQ/dr < 0, is rather thin. Given that the 
aspect ratio of magnetic flux belts is usually of order unity, it is difficult to envisage 
how one can explain the rather broad latitudinal distribution of flux of the same sign. 

As a possible solution to this problem one might think of the effects of magnetic 
helicity transport and possibly the anisotropy of the turbulent magnetic diffusivity. 
However, magnetic helicity fluxes modify primarily the effective a and it is not clear 
that they affect the aspect ratio of the toroidal flux belts. In the following we address 
the effect of an anisotropic magnetic diffusivity. Based on the dispersion relation for 
a dynamo wave in a two-dimensional domain (see, e.g., Brandenburg & Subramanian 
2005a, Sect. 6.5.2), 

A = - VT (k 2 x + ek 2 z ) + yJlaSk x , (23) 

where e is the degree of anisotropy, we find that in the marginally excited case, A = 0, 
the value of k x can be obtained iteratively. For this purpose, let us define the aspect 
ratio k x /k z , which is then given by 

k 2 + e = y/laSK/r&kl (24) 
We define the dynamo number as 

D = \aS/r&kl, (25) 
so equation ff24l reduces to 

n 2 = - e, (26) 

which can be solved iteratively. The result is shown in Fig. [15j In the special case e = 
we have n = D 1 ' 3 . 

It turns out that anisotropic diffusion tends to increase the ratio k = k x /k z , so the 
wavelength in the latitudinal (or x) direction becomes smaller, making the problem even 
worse. Obviously one must be careful with linear theory, so this result may not be too 
meaningful, but in the absence of any other evidence there is currently no particular 
reason to expect anisotropic diffusion being a solution to the problem of the aspect ratio. 

6.4- Predictability in distributed dynamos 

Flux transport dynamos have been used for predicting the strength of the next solar 
cycle, Cycle 24 (see Clarke 2006 for a general assessment). The outcome depends 
essentially on the way the observed solar activity is used to keep the evolution of the 
model in sync with the Sun. In Dikpati & Gilman (2005, 2006) the sunspot number is 
used as a proxy of the poloidal field production by a source term near the surface. On 
the other hand, in Choudhuri et al. (2007) the observed poloidal field around the poles is 
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Figure 15. Aspect ratio as a function of the degree of anisotropy, e for D — 2 (dashed 
line), 5 (dotted line), 10, 20, and 50. When the lines stop for small values of K = k x /k z 
there are no solutions any more. 




Figure 16. Simple dynamo model with dynamo number D — 4 with initial condition 
A = 0.01 and B — 0. The model reaches saturation at t/r ~ 3. The poloidal field is 
rescaled by factors 1.25, 1.23, 1.00, and 0.60 during the minima marked by vertical lines. 
The unsealed model is overplotted as a dotted line. Note the weak cycle amplitude at 
t/r w 5, corresponding to Cycle 24. 



used to correct the poloidal field in the upper part of their model during solar minimum. 
Regardless of the fact that the outcomes can be very different (strong Cycle 24 in the 
former model and weak Cycle 24 in the latter), it is clear that this topic has attracted 
significant attention (Clarke 2006, Tobias et al. 2006, Cameron & Schiissler 2006) and 
has led to the impression that the solar dynamo problem might be solved. 

The point of this section is to show that the predictive power of a model is not 
very sensitive to the details, and that even in the extreme case of a fully distributed toy 
model similar predictability can be achieved. To demonstrate this, we consider here the 
following simple model equations that are obtained from a single mode truncation of a 
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one-dimensional periodic model, i.e. 




(27) 



dB 
~dt 



ikVL'A - r~ x B 



(28) 



for the complex variables A and B, that characterize the poloidal and toroidal fields, 
respectively (Durney & Robinson 1982). Here, r = (rjtk)' 1 is the turbulent magnetic 
diffusion time, Q' is the radial gradient of the angular velocity, and a is the a effect. 
We assume here a simple form of a quenching with a = a /(l + \B\ 2 /B 2 q ), where a 
is the kinematic value of the a effect, and B eq is the equipartition field strength. The 
non-dimensional dynamo number, D = aofl' kLr~ 2 has to exceed the value D = 2 for 
dynamo action. 

Following Choudhuri et al. (2007) we use normalized values of the observed dipole 
moment (as given by Svalgaard et al. 2005) to correct the poloidal field amplitude A 
by factors 1.25, 1.23, 1.00, and 0.60 after Cycles 20, 21, 22, and 23, respectively. The 
result is shown in Fig. [16j The maxima in \B\ after each of the four rescalings are 1.17, 
1.17, 1.02, and 0.71 times the usual values. Applied to the Sun, this means that the 
next cycle will indeed be about 30% weaker than the previous one. This is obviously in 
perfect agreement with the model of Choudhuri et al. (2007) and of course the earlier 
investigations by Svalgaard et al. (2005). 

7. Conclusions 

There has been tremendous progress in dynamo theory since the beginning of the 
millennium, both theoretically and experimentally. In this paper we do not intend 
to review any of this progress. Instead, we have focussed on a number of new ideas that 
emerged in an attempt to connect theory and experiments. The relevance of experiments 
is not immediately obvious because in experiments the value of R m is currently below 
50 or so, while simulations can reach values on the order of 1000. On the other hand, 
in experiments the magnetic Prandtl number is small (10 -5 ), which is quite similar to 
the value in the Sun and other stars. To reach interesting values of R m that allow for 
dynamo action or other effects to occur, the fluid Reynolds number has to be very large 
(in excess of 10 6 or so), which is out of reach of simulations even in the intermediate 
future. 

On the other hand, experiments are now beginning to produce dynamo action in 
unconstrained flows (Monchaux et al. 2007), and corresponding experiments are now 
beginning to produce results that can meaningfully be compared with experiments. 
Examples include the work of Nore et al. (2006) for simulating flows close to those in 
the VKS experiment, and approximations to the VKS flow in terms of Taylor-Green 
flows (Ponty et al. 2004, Mininni et al. 2005b). Although neither of these cases is 
anywhere near the conditions relevant to the Sun, it is conceivable that the theoretical 
investigations that have been undertaken in just the last few years (e.g. regarding 
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magnetic helicity conservation and helicity fluxes) can lead to new paradigms that 
can in principle be tested experimentally (e.g. the phenomenon of inverse cascade-type 
behavior) . 

Even though progress is rapid, and success appears sometimes in reach, there 
are still a number of problems that are poorly understood. An example concerns the 
magnetic helicity flux whose effects are manifested in coronal mass ejections (Demoulin 
et al. 2002), but they have hardly been incorporated in dynamo models with the aim to 
alleviate the otherwise catastrophic quenching. 
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Appendix A. Magnetic versus kinetic helicity conservation 

In this appendix we want to show why there is such a dramatic difference between the 
conservation of magnetic and kinetic helicities. The kinetic helicity is indeed conserved 
if there is no magnetic field, i.e. no Lorentz force, and if v — exactly. Let us also 
assume V • A = V ■ U = for simplicity, although this is not critical. The pair of 
analogous equations can then be written as 
8 A 

— = UxB- v J-V<p + £ cxt , (A.l) 

?LL = Ux W-vQ-Vp + JF cxt . (A.2) 

For later reference we also quote here the curled evolution equations, where B = V x A 
is the magnetic field and W = V x U is the vorticity, i.e. 
dB 

= V x (U x B) -77V x J + V x £ ext , (A.3) 
= V x (U x W) - vV x Q + V x JF cxt . (A.4) 



dt 

dW 



dt 

Note that V-S = V- VF = 0is automatically fulfilled. Denoting, as usual, volume 
averages by angular brackets, and assuming periodic boundaries, the evolution equations 
for the two helicities are 

|-(A • B) = -2i](J ■ B) + 2(B ■ £ ext ), (A.5) 
dt 

lL(U ■ W) = -2u(Q ■ W) + 2(W ■ JFext). (A.6) 
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For completeness we give here also the evolution equations for the energy norms, 

|±(A 2 > = -(U ■ (A x B)) - V (B 2 ) + (A ■ S cxt ), (A.7) 



|±(C/ 2 > = - (17 (UxW)) -v(W 2 ) + (U ■ T, xl ). (A.8) 



^\{B 2 ) = -(17 • (J x B)> - ^(J 2 ) + (J • £ ex t), (A.9) 



j t \(W 2 ) = -(U • (Q x W)> - z/(Q 2 > + (Q ■ JP cxt ). (A.10) 



i-(A • B) = - (U (BxB)) - r)(J -B) + (B- £ cxt ) (A.ll) 
2 X J7 . W) = - (U (WxW)) - u{Q ■ W) + (W ■ JF cxt > (A.12) 



Note that we kept the term {U ■ (U x W)) = in order to enhance the formal analogy 
of the equations. In the helicity equations one could restore two similarly redundant 
terms, 

d_ 

dt 

In conclusion, all the evolution equations have the following three terms: an internal 
driving term such as (U • (J x B)), an external driving term such as (U ■ tF ext ), and 
a viscous or resistive loss term. In the magnetic energy equation there is the internal 
driving term (work done against the Lorentz force), but usually no external driving 
term, since £ cxt = is assumed here. The kinetic energy equation, on the other hand, 
does not have an internal driving term, because (U ■ (U x W)) = 0. Thus, both kinetic 
and magnetic energy equations have driving terms, so we have 



|i(B 2 ) = -(17- (J xB))- v (JZ), (A.13) 



|i(t/ 2 ) = (C/.^ cxt ) - u(W'). (A-ll) 



Thus, in a statistically steady state, the dissipative terms have to balance the 
corresponding driving terms. In the limit of large magnetic and fluid Reynolds numbers 
this leads to the asymptotic scalings 

|J|~?T 1/2 , \W\^iy- l/2 . (A.15) 

At the same time, because kinetic and magnetic energies are bounded, magnetic field 
strength and velocity do not diverge but stay independent of r\ and u, respectively. This 
also implies that the typical inverse length scales scale like 

k ~ |W|/|J7| ~ z/~ 1/2 , (A.16) 

and hence that 

\Q\ „ k\W\ ~ u~\ (A.17) 
so that the magnetic and kinetic helicity dissipation terms scale like 

\i](J ■ B)\ — > i] 1 / 2 — > 0, (A.18) 
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\u(Q- W)\ -> u- 1/2 -> oo, (A.19) 

which highlights the fundamental difference between kinetic and magnetic helicity 
conservation when 77 and v are not zero but small. 

Appendix B. The forcing function 

For completeness we specify here the forcing function used in the present papejjj. It is 
defined as 

f(x, t) = Re{Nf k{t) exp[ifc(t) • x + i0(t)]}, (B.l) 

where x is the position vector. The wavevector k(t) and the random phase — tt < <p(t) < 
7r change at every time step, so f(x,t) is (^-correlated in time. For the time-integrated 
forcing function to be independent of the length of the time step St, the normalization 
factor N has to be proportional to St~ 1 ^ 2 . On dimensional grounds it is chosen to be 
N = fQC s (\k\c & / St) 1 / 2 , where /o is a nondimensional forcing amplitude. At each timestep 
we select randomly one of many possible wavevectors in a certain range around a given 
forcing wavenumber. The average wavenumber is referred to as kf. We force the system 
with transverse helical waves, 



t ohel) with R = Jtnvk h 
where a = 1 for positive helicity of the forcing function 



f k = R -fk with ^ = Vi+X ' (R2) 



/^ nohel) = (k x e) /Jk*-{k-e)\ (B.3) 

is a non-helical forcing function, and e is an arbitrary unit vector not aligned with k; 
note that |/^,| 2 = 1. 

Appendix C. One-mode truncation of the Blackman— Field model 

A one-mode truncation of the dynamically quenched a effect has been studied by 
Blackman & Brandenburg (2002), where B(z,t) = B(t) exp(ikxz) has been assumed. 
Here we include the explicit time dependence of the electromotive force £(z,t) = 
£(t) exp(ikiz), so our model equations are 

^L = ik 1 x(£-r ! J), (C.l) 
at 

^- = aB-fj t J--, (C.2) 



dt ' t 



da 



dT - - 2 «*> 



Re(C'-B) q M 

- D eq - rL m 



(C.3) 



|| This forcing function was also used by Brandenburg (2001), but in his Eq. (5) the factor 2 in the 
denominator should have been replaced by y/2 for a proper normalization. 
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where B and £ are complex and «m a real dependent variable, k\ = (0, 0, k\) T is the 
wave vector, a = % + ctyi with a = ret and r/ t = T^ t , and r is the relaxation time. The 
closure assumption consists in representing the triple correlations by the damping term 
S/r. 

In order to associate this model with a simulation, we use the values of kf/ki, 
R m and chose a value of St to determine the following set of model parameters: 
M rms = 3r/ t fcf/St, r = St/(-u rms /cf), and ax = Vth, which assumes perfectly helical 
turbulence, which is obviously not realistic for Taylor- Green flow dynamos. 



